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Abstract 

The pion nucleon vertex function at finite temperature is studied in the framework 
of: (a) the thermal (linear) sigma model to leading (one-loop) order, and (b) a 
thermal QCD-Finite Energy Sum Rule. Results from both methods indicate that 
the strength of the pion-nucleon coupling decreases with increasing T, vanishing 
at a critical temperature. The associated mean-square radius is a monotonically 
increasing function of T, diverging at the critical temperature. This is interpreted as 
(analytical) evidence for the quark-gluon deconfinement phase transition. 



The temperature behaviour of hadronic Green's functions, and their associated parameters 
such as masses, widths, couplings, etc., has received considerable attention lately, given its 
impact on the search for the quark-gluon plasma Two successful theoretical frameworks 
for these studies are the thermal sigma model |2|]-[|], and QCD sum rules ||-||. The 
former technique provides information on the T-dependence of pion and nucleon masses 
and widths associated, respectively, with the real and imaginary parts of their two-point 
Green's functions. While these masses show no appreciable variation with temperature, 
||-[[| their widths exhibit a dramatic increase with increasing T |],||. This result is in 
line with the expectation that hadronic widths, interpreted as absorption coefficients in 
the thermal bath, should diverge at some critical temperature |7J. This provides a signal 
or phenomenological order parameter for the quark-gluon deconfmement phase transition. 
Another such signal is the thermal behaviour of hadronic couplings and form factors (three- 
point functions), which should vanish at a critical temperature, where the associated mean 
square radii should diverge. This has been explicitly confirmed for the electromagnetic 
form factor of the pion 0, and for the rho-pi-pi coupling 11. In this note we study the 
irNN vertex function at finite temperature using the thermal (linear) sigma model, as well 
as QCD sum rules. The purpose is to obtain additional confirming (analytical) evidence for 
the deconfmement phase transition, as well as information on this vertex function, which 
should be of use in hadron gas models at finite temperature. 

We begin with the (linear) sigma model, and consider the 7rNN vertex 

IV) = V(q 2 )u f (p') l5 T aUl (p) (1) 
where the nucleons are on-mass shell, and the pion has virtual mass q 2 = (p' — p) 2 . The 



renormalization of the sigma model (at T = 0) is discussed e.g. in [ID , and the renormal- 
ization of the 7rNN vertex (before the invention of dimensional regularization and the MS 
scheme) may be found in [|ll|]. At finite temperature we shall use the Dolan-Jackiw real 



time propagators ||12|| , together with the fact that thermal corrections do not induce any 
new kind of ultraviolet corrections. Hence, the thermal theory may be renormalized as at 
T = 0. We have done this using dimensional regularization and the MS scheme. To lead- 
ing (one-loop) order, the relevant diagrams are shown in Fig. 1 (a-e). At the kinematical 
point q 2 =0, and in the chiral limit, the expression for the irreducible vertex V(q 2 ) is given 



by 



where 
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with Mjy and M a being the nucleon and sigma meson masses, respectively. Equation (2) 
may be regarded as an expression for the effective coupling constant g^NN i n the chiral 
limit. It is equal to g if /3(0) vanishes. This happens for M a ~ 1300 MeV; in fact, /3(0) is 
small and negative if M a is bigger than this value. One should recall that in the linear sigma 
model gA = 1, so that the Goldberger-Treiman relation (GTR) |0J becomes: M^ = gf n . 
Using the chiral symmetry limit values fl4| : Mjv — 800 MeV, and f n ~ 80 MeV, one 
finds: g ~ 10, not far from the experimental value g^NN — 13. In any case, here we are 
only interested in the temperature behaviour of the ratio V(q 2 , T) /V(q 2 , 0); particularly in 
the possibility that this ratio vanishes at a critical temperature, and that the mean square 
radius diverges there. This will turn out to be largely independent of the particular value 
assumed by g or g n NN and /5(0) (or equivalently M a ), although the specific value of the 
critical temperature does depend on the latter. 



Turning to the temperature corrections to the graphs shown in Fig. 1, we need only consider 
Fig. 1(b) and Fig. 1(d) (identical to Fig. 1(e)), as Fig. 1(c) is Boltzmann suppressed on 
account of Mjy — M a » fi^. The thermal correction to the graph Fig. 1(b) is given by 



w 2 rps 3- / a ( \ [ dAk k 2 n B (k )S(k 2 - /4) 

r(g ,T) = -g u f ( P ) W,(P) J W [{pf _ W _ m [(p _ W _ ^ (4) 

where ub is the Bose thermal factor: ns(z) = [e z l T — l) -1 . Hence, this contribution 
vanishes in the chiral limit. That of Fig. 1(d) is found to be 
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We choose for convenience a Lorentz frame in which the incoming nucleon is at rest with 
respect to the heat bath ((p = 0). We have checked that the final results are largely 
independent of the choice of frame. The thermal effective irNN coupling in the chiral 
limit, and at q 2 =0, may then be written as 



mr) = 1 _ 9 2 t 2 

V(0,0) 12M 2 (l + <? 2 /3(0)) U 

where the GTR has been used. Notice that an extrapolation in temperature of this result 
implies a critical temperature 



T d = ^12(l + (3(0)9 2 )M N /g (7) 



d ^ 



which depends on the value of the sigma-meson mass through (3(0); numerically, T ( 
150 — 300 MeV if M a ~ 1300 — 1600 MeV . One should not assign too much importance to 
the specific numerical values of this critical temperature and/or the sigma-meson mass, to 
wit. First, the relations among parameters in the sigma model are valid at the 25-30% level; 
e.g. qa = 1 instead of the experimental value Qa = 1-26, and q-^nn — 10 (from the GTR), as 
opposed to the experimental value g n NN — 13, etc.. Second, Eq.(6) and the T 2 dependence 
is a consequence of the one-loop approximation. Higher loop corrections will induce higher 
order (in T) terms which will alter the numerical value of the critical temperature. These 
will be suppressed, though, by inverse powers of the nucleon/sigma-meson masses. A 
similar situation arises in chiral perturbation theory and the T-dependence of the pion 
decay constant / T (T). To order T 2 , f n (T) vanishes at T c ~ 240 MeV fL5f , while higher 



order corrections bring down this value considerably. What we find important here, is that 
the 7iNN coupling at leading order in T decreases with increasing temperature. 

Next, we consider the mean-square radius associated with V(q 2 ,T), and defined as 
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<rl NN > T _ V(0,0) 
<rl NN>0 1/(0, T) 
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dq 2 



dq 2 



(8) 



We have calculated this ratio numerically, the result being shown in Figure 2. An extrap- 
olation in temperature indicates quite clearly the divergence of this mean square radius. 
This may be interpreted as a signal for deconfinement, to the extent that the size of the 
nucleon, as probed by a pion, increases with increasing temperature, becoming infinite at 
T = T d . 



We study next the same vertex function, but in the framework of QCD Finite Energy 
Sum Rules (FESR). This will provide important independent support to the above result, 
especially since the QCD sum rule technique, unlike the sigma model, does not entail 
any expansion in powers of the temperature. To this end, we begin the analysis at zero 
temperature and introduce the three-point function 



U(p,p',q) =i 2 J J d 4 x d A y{0\T(r](x)J 5 (y)f}(0))\0)e i(p '- x 



-Q-y) 



(9) 



where the nucleon and pion interpolating currents, r)(x) and J*,(x), respectively, are chosen 

as 

(10) 
(11) 



•nix) = e abc [u a {x)C^u b (x)] ^ 5 d c (x) 
■h(x) = i u(x)j5u(x) — d(x)j5d(x) 



where C is the charge conjugation operator. The couplings of these currents to the nucleon 
and the pion are defined as 



(0\r ] (0)\N(p,s)) = \ N u(p,s) 
{N(p 2 , s 2 )| J 5 (0)|7V(pi, si)) = u(p 2 , s 2 )l59p{q 2 )u{pi, si) 



where 



2 . fir^l 9-kNN 



m q q 2 - /4 



(12) 
(13) 

(14) 
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and where m q is the average of the up and down quark masses. In our normalization, 
the pion decay constant is f n ~ 93 MeV. The hadronic representation of the imaginary 
part of the vertex function Eq. (9) is obtained by inserting a complete set of hadronic 
states. After summing over spins, and making the standard nucleon-pole approximation 
(thus including in the continuum all the radial excitations of the nucleon) one obtains 



Since we are interested in the pion-nucleon coupling in the vicinity of q 2 =0, we can safely 
neglect any q 2 dependence in g-^NN- This dependence would arise from the contribution of 
the radial excitations of the pion, 7r'(1300) etc., which in the chiral limit is a correction of 
order q 2 /M 2 ,. As usual, the hadronic continuum with thresholds s and s' is modelled by 
the QCD spectral function. The leading order diagrams needed to compute the latter are 
shown in Fig. 3. In the chiral limit, the relevant structure to be sought is proportional to 
qlq 2 . It turns out that the diagram Fig. 3a does not have this behaviour, while that of Fig. 
3b (plus all other related diagrams) gives 



where < uu >~< dd >=< qq > has been used. By means of Cauchy theorem, and 
assuming global (quark-hadron) duality, one obtains the lowest dimensional FESR 



ImU(s,s',q 2 )\ EAD 



g P X 2 N M N (* 75 4) n 2 5(s - M 2 N )5(s' - M 2 N ) 
+ e( s - So )0( s ' - s' )ImU(s, s', g 2 ) |qcd 



(15) 




(16) 



9ttNN — 
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(17) 



where use has been made of the Gell-Mann, Oakes and Renner relation JRI] 




2m g (qq) 



(18) 
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Since the dispersion in p 2 = s and p' 2 = s' refers to the two nucleonic legs, it is reasonable 
to assume sq — s'q. An analysis of the two-point function involving the nucleonic current 
r](x) |I7| in the framework of QCD FESR yields 



An ~ 192tt 4 



X 2 N M N 



(qq) „ a 

r 2 6 



8tt 2 



(19) 



which determines the nucleon mass in terms of sq. Conversely, using Mat and < qq > as 
input, Eq.(19) fixes and s . In this fashion, Eq.(17) becomes: g n NN = 487r/ 7r /Mjv — 15, 
not far from the experimental value g n NN — 13. This level of agreement is more than 
enough for our purpose here, which is to determine the temperature behaviour of the 
tiNN coupling, i.e. the ratio (^ata^T) #^^(0). To achieve this, we have calculated the 
thermal corrections to the QCD spectral function with the result 



ImU(p,p',q) = ^(m)-^r {p 2 f(P>T)+p' 2 f(p',T)} 



where 



(20) 



f(p,T) 



dx 



1 — n F 
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(21) 



and a similar expression for f(p',T). The FESR in this case becomes 



where use has been made of the thermal Gell-Mann, Oakes and Renner relation, which has 
been recently shown to be valid over a wide range of temperatures |[L8|| . The temperature 



behaviour of /„., valid up to the critical temperature, has been obtained in JO]]. The 
function sq(T) has been determined from a FESR for the two-point function involving the 
axial-vector current p3 - [[H]]; it scales as: so(T)/s (0) — (/ 7r (T)// 7r (0)) 2 . The thermal 
nucleonic coupling Xn(T) has been determined from a QCD-FESR in the nucleon channel, 
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with the result Wl 



A* CO = A^(0) 



(23) 



where 



57b r\/ s o(-L) ru/* ru/-z 
G a {T) = - — ——3 / dw dx dy x(u - 2x) {-n F (x) 

(S (T)) J J Jlu/2-x 

-n F (y) + n F (x)n F (y) + n F (^ - x - y) [n F {x] + nj?(j/) - 1]} . (24) 



and n F (z) = (1 + e - */ 7 ) -1 is the Fermi factor. With all this information one can then 
solve Eq.(22), after choosing a particular Lorentz frame. Our choice is the rest frame of 
the incoming nucleon, i.e. p = and po = y/s, however, the final results are quite insen- 
sitive to the choice of frame. The result for the ratio g n NN(T) / g nNN (0) as a function of 
T/Trf is shown in Fig. 4. One can clearly appreciate the vanishing of this coupling at the 
critical temperature. This is the temperature at which f w (T) vanishes, i.e. the critical 
temperature for the chiral-symmetry restoration phase transition, which is basically the 
same temperature at which Sq(T) vanishes, i.e. the critical temperature for the quark-gluon 
deconfinement phase transition P0| - [2"T|. 



Finally, the mean-square radius associated to the pion-nucleon vertex 



d 

t Inn)t = 6— \ng n NN(q 2 ,T)\ q 2 =0 



dq 



(25) 



can be easily calculated, with the result 



r s (T) rs' {T) rl 

J J ds ds' J dx [s (l - 2n F (v^/2)J + s' (1 - 2n F (z)) 



-1 



x— I f° dsds r -^= [ dxn 2 F (z)e z / T 
2T jo Jo vs 7-i 



1 + x 



s + s' 



\s — s 



where 



s + s' + \s — s'\x 



(26) 



(27) 



and the rest frame of the incoming nucleon has been used. A numerical evaluation of 
Eq.(26) gives the result shown in Fig. 5. Notice that since we have made the pion-pole ap- 
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proximation, g^NN at T = is independent of q 2 . The mean-square radius is non- vanishing 
only at finite temperature, where a q 2 dependence appears through the Fermi factors. 

In summary, the vanishing of the pion-nucleon coupling, and the divergence of the associ- 
ated mean-square radius, at a critical temperature has been shown to follow from the ther- 
mal linear sigma model at leading (one-loop) order, as well as from a thermal QCD-FESR. 
This may be viewed as (analytical) evidence supporting the existence of the quark-gluon 
deconfinement phase transition. As the critical temperature is approached, the strength of 
the coupling of pions to nucleons is quenched, and at the same time, the size of the nucleon 
as probed by the pion gets bigger. The qualitative agreement between the two methods 
lends further support to the extension of the QCD sum rule program to finite temperature. 
It should be noticed that potential non-diagonal vacuum condensates |22| do not enter our 
FESR because of their (higher) dimensionality. We have emphasized many times in the 
past that QCD-FESR are far better than e.g. QCD-Laplace transform sum rules, to the 
extent that the lowest dimensional thermal FESR do not involve (unknown) non-diagonal 
vacuuum condensates. 
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Figure Captions 



Figure 1. Leading (one- loop) diagrams contributing to g n NN in the linear sigma model. 
Figure 2. Thermal behaviour of the irNN mean square radius, Eq.(8), in the linear sigma 
model. 

Figure 3. Leading order QCD diagrams entering the determination of the spectral function 
relevant to g n NN- 

Figure 4. Thermal behaviour of the 7rNN coupling determined from the QCD-FESR 
Eq.(22). 

Figure 5. Thermal behaviour of the ttNN mean square radius, Eq.(26), according to the 
QCD-FESR. 
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